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Abstract 

Magnetic Heusler materials with very low Gilbert damping are expected to show novel magnonic 
transport phenomena. We report nonlinear generation of higher harmonics leading to the emission 
of caustic spin-wave beams in a low-damping, micro-structured C02Mno.6Feo.4Si Heusler waveguide. 
The source for the higher harmonic generation is a localized edge mode formed by the strongly 
inhomogeneous field distribution at the edges of the spin-wave waveguide. The radiation character- 
istics of the propagating caustic waves observed at twice and three times the excitation frequency 
are described by an analytical calculation based on the anisotropic dispersion of spin waves in a 
magnetic thin film. 
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In the last years nonlinear spin dynamics in magnetic micro-structures made of metallic 
ferromagnetic thin films or layer stacks have gained large interest. l|-[5( The intrinsically 
nonlinear Landau-Lifshitz and Gilbert equation (LLG), that governs the spin dynamics, 
gives rise to a variety of nonlinear effects, (g, |7] 

Among the metallic ferromagnets, the class of Cobalt-based Heusler materials is very 
promising for future magnon spintronics devices and the observation of new phenomena of 
magnonic transport. The reasons for the interest in these materials are the small magnetic 
Gilbert damping, the high spin-polarization, and the high Curie temperature. [8] 

As shown recently, the full Heusler compound C02Mno.6Feo.4Si (CMFS) is a very suitable 
material to be used as a micro-structured spin-wave waveguide due to the increased decay 
length which was observed for wave propagation in the linear regime. {9] The reason for this 
observation is the low Gilbert damping of a = 3 x 10 -3 of CMFS compared to Ni 81 Feig with 
a = 8 x 10" 3 , which is the material commonly used in related studies. js] The decreased 
magnetic losses not only lead to an increase of the decay length but also to large precession 
angles of the magnetic moments and, thus, to the occurrence of nonlinear effects. Regarding 
future applications, the investigation and the thorough understanding of phenomena related 
to the spin-wave propagation in the nonlinear regime in Heusler compounds is crucial. 



In this Letter, we report nonlinear higher harmonic generation from a localized edge 
mode causing the emission of caustic spin-wave beams [12| , 



in a micro-structured 



Heusler waveguide. Spin-wave caustics are characterized by the small transversal aperture 
of a beam, which practically does not increase during propagation, and the well-defined 
direction of propagation. 

The investigated sample is a 5 /im wide spin- wave waveguide structured from a 30 nm 
thick film of the Heusler compound CMFS. Details about the fabrication and the material 



properties can be found in Refs. 



14 



and 
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The microfabrication of the waveguide was per- 
formed using electron-beam lithography and ion-milling. For the excitation of spin dynamics 
in the waveguide the shortened end of a coplanar waveguide made of copper was placed on 
top of it. The Oersted field created by a microwave current in this antenna structure can be 
used to excite spin dynamics in the Gigahertz range. The antenna has a thickness of 400 nm 
and a width of Ax = 1 /im. 

All observations have been carried out using Brillouin light scattering microscopy 
(/iBLS).[16] /iBLS is a powerful tool to investigate spin dynamics in microstructures with a 



antenna structure and 
exciting field 




external 
field 



FIG. 1. (Color online) Sketch of the sample design. The shortened end of a coplanar waveguide 
is used as an antenna structure to excite spin dynamics in a 5 /xm wide CMFS waveguide with a 
thickness of 30 nm. The waveguide is positioned in the x-y-plane with its long axis pointing along 
the x-direction. The external magnetic field is applied transversely to the waveguide in y-direction. 
The figure includes a /iBLS spectrum taken at a distance of 4.5 //m from the antenna in the center 
of the waveguide (see laser beam in the sketch) for an excitation frequency of f e = 3.5 GHz, a 
microwave power of 20 mW, and an external field fioH ex t = 48 mT. 



spatial resolution of about 250 nm and a frequency resolution of up to 50 MHz. 

In the following description, the waveguide is positioned in the x-y-plane with the long 
axis pointing in x-direction. The origin of the coordinate system is given by the position 
of the antenna between x = — 1 /im and x = /im. An external magnetic field of /j,oH ext = 
48 mT was japplied transversely to the waveguide in y-direction resulting in Damon-Eshbach 
geometry |17| for spin waves propagating along the waveguide. A sketch of the sample layout 
is shown in Fig.[TJ 

In addition, Fig.[T] includes a spectrum taken by /xBLS for an excitation frequency of 
f e = 3.5 GHz and a microwave power of 20 mW at a distance of 4.5 /im from the antenna in 
the center of the waveguide. The spectrum shows not only a peak at f e = 3.5 GHz but also 
at 2/ e = 7.0 GHz and 3/ e = 10.5 GHz. Furthermore, the intensity of the directly excited spin 
wave at 3.5 GHz is lower than for the higher harmonics at the point of observation. As we will 
see, the higher harmonics are excited resonantly by nonlinear magnon-magnon interactions 
and the intensity distribution is a consequence of the different propagation characteristics 
of the observed spin-wave modes. 



A two-dimensional intensity distribution 
as well as the calculated dispersion relation 



or the detection frequency fd = fe = 3.5 GHz 



18] for the center of the CMFS waveguide are 
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FIG. 2. (Color online) (a) /^BLS intensity distribution for the detection frequency /a being equal 
to the excitation frequency f e = 3.5 GHz. The observed intensity of the edge mode is maximal near 
the edges of the CMFS waveguide. Please note that the exciting antenna is positioned between 
x = — 1 /im and x = fim. The circle in the graph indicates the position of the measurement 
presented in Fig.Jjl (b) Calculated dispersion relation for the center of the CMFS waveguide 
according to Ref. 18 as well as excitation frequency and higher harmonics (dashed lines). 



shown in Fig.[2j Figure[2](a) reveals a strong localization of the intensity at the edges of the 
waveguide. The non- vanishing intensity close to the antenna and between the edges of the 
waveguide (y = 1 — 4/im) can be attributed to nonresonant, forced excitation by the Oersted 



field created by the microwave current. Figure[2[ 



d) shows the spin-wave dispersion for the 



CMFS waveguide calculated according to Ref.[18| assuming a homogeneous magnetization 
oriented in y-direction by the external field. As can be seen, the lower cut-off frequency 
/ = 6.9 GHz is well above the excitation frequency of f e = 3.5 GHz (see dashed line in 
Fig.[2](b)). In the center of the waveguide, the assumption of a homogeneous magnetization is 
a very good approximation. However, the demagnetizing fields created by the magnetization 
pointing transversely to the waveguide are responsible for a strongly inhomogeneous and 
decreased effective field near its edges. This field configuration allows for the existence of 
localized spin waves - commonly referred to as edge modes - energetically far below the spin- 
wave dispersion for propagating modes in the waveguide, 1CJ, 111) Therefore, we conclude that 
the spin-wave mode at fd = fc = 3.5 GHz is excited resonantly by the microwave field. 

As can be seen from the spectrum in Fig.fi], the applied microwave power of 20 mW is 
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FIG. 3. (Color online) Power dependence of the BLS intensity for the detection frequencies fd = 
3, 5 GHz, 7 GHz, and 10.5 GHz for the fixed excitation frequency f e = 3.5 GHz. Please note the 
log-log presentation of the data. The lines in the graph correspond to fits according to Eq.[TJ A 
least square fit of the data yields sif = 0.9 ± 0.1, S2f = 2.1 ± 0.1, and S3f = 2.8 ± 0.3. 

sufficiently high to observe the nonlinear generation of higher harmonics of the excitation 
frequency f e in the /zBLS spectra. These resonant frequency multiplications to / = 2f c = 
7.0 GHz and / = 3f e = 10.5 GHz result in the excitation of propagating spin-wave modes in 
the waveguide energetically above the cut-off frequency of the dispersion shown in Fig.^b). 

Figure|3] shows the dependence of the directly excited mode and the higher harmonics 
on the applied microwave power. This data has been acquired close to the position of the 
edge mode and near the antenna at x = 0.7 /im and y = 0.8 /im. The data is presented on 
a log-log scale with fits according to 



I n (p) = A nP s - + b, (1) 

where I n is the /xBLS intensity, A n a coupling parameter, p the applied microwave power, 
and b the noise- level in our measurement. As expected, these processes do not show a 
threshold power level, but reliable detection on the background of the noise is not possible 
for powers below 5mW for = 3/ = 10.5 GHz. The different slopes of the curves for the 
different spin-wave modes n are caused by the different power-laws specified by the exponent 
s n . A least square fit of the data yields sn = 0.9 ± 0.1, S2f = 2.1 ± 0.1, and S3f = 2.8 ± 0.3. 
These experimental findings close to the integer values 1, 2, and 3 are in accordance with 
both reported experimental data and theoretical predictions for the nonlinear generation of 



higher harmonics. 

The observation of the second harmonic can be understood qualitatively by considering 
the strong demagnetizing fields caused by the out-of-plane component m z (t) during the 
magnetization precession. Due to the demagnetizing fields, the magnetization precession 
M(i) around the y-direction (defined by the effective field) follows an elliptical trajectory 
rather than a circular one. In contrast to the case of a circular precession, the resulting 
projection of M on the y-axis is time dependent and oscillating with the frequency 2/ e . The 
resulting dynamic dipolar field |h y (t)| oc m 2 x — m 2 z can be regarded as the source for the 
frequency doubling. Similar considerations lead to the observation of higher harmonics. A 
full quantitative derivation of higher harmonic generation and other nonlinear effects based 
on the expansion of the LLG in terms of the dynamic magnetization can be found in Ref.Q. 

FiguresHJa) and (b) show intensity maps for the detection frequencies /a = 2f c = 7.0 GHz 
and /d = 3/ e = 10.5 GHz, respectively. In both cases the intensity radiated from the 
position of the edge mode is strongly directed, has a small transversal aperture, and shows 
nondiffr active behavior. These radiation characteristics recorded for the three spin-wave 
modes at / = 3.5, 7.0 and 10.5 GHz presented in Figs.|2^a) and H] are responsible for the 
intensity distribution shown in the spectrum in Fig.[TJ The position of this measurement 
is indicated in the corresponding intensity maps with a circle. Since the spin wave at 
/ e = 3.5 GHz is localized at the edges of the waveguide, its intensity is comparably weak in 
the center. In contrast, the higher harmonics have frequencies above the cut-off frequency of 
the spin-wave dispersion and can propagate in the center of the waveguide. The propagation 
angles of these modes support the intensity distribution recorded in our measurement. While 
an increased intensity can be found already for /d = 2/ e , for /d = 3/ e the two beams starting 
from both edges of the spin-wave waveguide even intersect at the measurement position 
resulting in the highest intensity at this point. Because of the well-defined propagation angles 
of the higher harmonics and the localization of the edge mode, the intensity distribution is 
strongly depending on the measurement position. 

This observation of spin-wave beams with small transversal aperture is reminiscent of the 



results in Refs. 



19 



and 



20 



where nonlinear three-magnon scattering in yttrium iron garnet 
is reported. However, in that case, the propagation direction of the nonlinearly generated 
spin wave is given by momentum conservation in the scattering process. In contrary, in 
our case, due to the strong localization of the edge mode, the assumption of a well-defined 




FIG. 4. (Color online) /zBLS intensity distribution for (a) fd = 2/ e = 7.0 GHz and (b) /d = 
3/ e = 10.5 GHz. Both intensity maps show strongly directed spin-wave beams along the angle 
9 = Z(H ext , vq). The lines in the maps are guides to the eye to identify the propagation angle 6. 
The circles in the graphs indicate the position of the measurement presented in Fig.[TJ 

initial wavevector and, thus, a strict momentum conservation is not justified. In particular, 
it is not possible to find an initial wavevector that allows for the nonlinear generation of the 
second and third harmonic at the same time still respecting momentum conservation. 

In the following, we will describe the observed propagation characteristics using the prop- 
erties of the anisotropic spin-wave dispersion in a magnetic thin film. jl2. 13 . 18 ] Because 
of this anisotropy, the direction of the flow of energy, which is given by the direction of 
the group velocity vg = 27r<9/(k)/<9k of the investigated spin waves, can differ significantly 
from the direction of its wavevector k. To estimate the relevant range in k-space in our 
experiment, we have to consider the lateral dimensions of the source for the nonlinear pro- 
cesses. Since the excitation by the oscillating Oersted field is most efficient directly below 
the antenna, the edge mode has the highest intensity in this region given by the width of 
the antenna of Ax = 1 fim. The spread of the edge mode in y-direction can be estimated 
from the intensity map in Fig.[2^a) to be smaller than 1 /im. Because of this localization, the 
edge mode in our measurement at / = 3.5 GHz can be regarded as a source for the nonlinear 
emission of the higher harmonics with lateral dimensions of approximately lxl /im 2 . As a 
first approximation, the Fourier transformation of this geometry lets us estimate the maxi- 
mum wavevector that can be excited by the edge mode to be k max 6.3 /im" 1 . As we will 
see, the direction of the group velocity can be assumed to be constant for most wavevectors 



that can be excited. This finally leads to the formation of the caustics in our experiment. 

For given frequency and external field, the iso-frequency curve f(k x , k y ) = const can be 
calculated analytically from the dispersion relation. Calculations for / = 2f e = 7.0 GHz and 
/ = 3/c = 10.5 GHz are illustrated in Fig.^a), where k y is shown as a function of k x . Using 
this data, we calculate the direction 9 of the flow of energy of the spin waves relative to the 
externally applied field by: 

9 = Z(H ex t, vg) = arctan(f 2/ /t> ;c ) = arctan(d£; y /<i&; x ). (2) 

Figure[5]^b) shows the calculated propagation angle 9 in the CMFS waveguide as a function 
of k y . The most important feature in the trend of 9 is the small variation of A9 < 2° in the 
range of k y = 2 — 7/im -1 for both frequencies / = 2/ e and / = 3/ e . While the wavevector 
k changes, the direction of vq - and, thus, the flow of energy - keeps almost constant as 
a function of k. In this range, which includes the maximum wavevector /c max ~ 6.3 /im -1 
that can be emitted from the edge mode (see considerations above), the calculations yield 
0caic(2/) = 79° and C aic(3/) = 66° as mean values, respectively. The dash-dotted lines in 
Fig.0(b) represent the propagation angles 9 of the spin-wave beams observed experimentally 
as shown in Fig.H] (9 exp (2f) = 78° and 9 exp (3f) = 67°). The comparison of experimental 
findings and analytical calculations shows an agreement within the expected accuracy of 
our measurement setup and is, therefore, supporting our conclusion. Higher harmonics with 
k y < 2 /im -1 are emitted with strongly varying directions from the edge mode and can be 
regarded as a negligible background in our measurement. 

In summary, we reported nonlinear higher harmonic generation from a localized source 
in a micro-structured CMFS waveguide leading to the emission of strongly directed spin- 
wave beams or caustics. This observation results from the complex interplay of different 
phenomena in magnonic transport in magnetic microstructures. The localization of an edge 
mode due to demagnetizing fields in the waveguide leads to the formation of a source for 
the following nonlinear processes. The nonlinear higher harmonic generation results in the 
resonant excitation and emission of propagating spin waves at 2/ e and 3/ e in a wavevector 
range corresponding to the localization of the edge mode. The experimentally observed 
power dependencies of the different spin-wave modes show the expected behavior for direct 
resonant excitation and nonlinear higher harmonic generation. As shown by our calculation, 
the anisotropic spin-wave dispersion yields a well-defined direction of the flow of energy 
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FIG. 5. (Color online) Analytical calculations according to Ref.Q- (a) Iso-frequency curves 
f const = /(t, k y ) for / = 2/ e = 7.0 GHz and / = 3/ c = 10.5 GHz. Based on these calculations ex- 
emplary directions for k, vq, H ex t and the propagation angle 9 = Z(H ex t, Vq) = arctan(d/c y /d/c x ) 
are shown in the graph, (b) Radiation direction 6 calculated from the iso-frequency curves shown 
in (a). Dash-dotted lines correspond to the angles 9 exp observed in the experiment. 



of the emitted spin waves in the relevant range in k-space. The calculation is not only 
qualitatively in accordance with our experimental findings but does also show quantitative 
agreement. 
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